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1 Introduction 



The direct approach to investigate integrable Hamiltonian systems is a very classical subject 02]. It consists 
in determining the class of potentials supporting additional invariants within some specified family of phase- 
space functions. This method produced several interesting results in the eighties, as illustrated in the review 
by Hietarinta j3] ■ Recently 0J 03 E] the approach has been applied to treat in a unified way both invariants 
at fixed and arbitrary energy. 

Many results are known for natural reversible Hamiltonians. One of the reasons for this is that the search 
for additional invariants can be restricted to functions with a definite parity in the momenta. This property 
leads to a substantial reduction in the usually very complicated set of equations. Much less is known in the 
case of Hamiltonians with vector potentials. For a long time, the only systematic attempt to cope with this 
case was that of Dorizzi, Grammaticos, Ramani and Winternitz [7] providing a set of solutions in Cartesian 
coordinates. Recently, McSween and Winternitz 8 obtained some new solutions in polar coordinates and 
Berube and Winternitz |!|| extended the results to the corresponding quantum problem. In both works the 
authors also identify the subset of superintegrable systems. In an attempt to extend these results to include 
weak integrability, we have provided 1 1 1 1 a general solution for linear invariants and analyzed some new 
classes of weakly-integrable systems. 

The purpose of the present paper is to reinvestigate Hamiltonian systems with both scalar and vector 
potentials, trying to identify those admitting the existence of a second invariant which is a quadratic poly- 
nomial in the momenta. We state the general approach at arbitrary and fixed value of the first invariant 
(Jacobi constant) and show that, in the case of strong integrability, it is possible to get a general formal 
solution. This is valid for every standard coordinate systems which are the same as the separable ones in 
the purely scalar case. Therefore, in addition to the above mentioned Cartesian and polar case, solutions 
in parabolic and elliptical coordinates can be looked for. In all cases, the potentials are defined in terms of 
a pair of scalar functions for which we get the integrability conditions: solving them determines the vector 
potential, whereas the scalar potential is subject to an additional linear differential equation. We provide 
some new examples of integrable systems with a vector potential whose existence can be discovered working 
in parabolic and elliptical coordinates. 

The plan of the paper is as follows: in section 2 we recall the structure of Hamiltonian systems with terms 
linear in the velocities; in section 3 we illustrate a version of the direct approach to find polynomial invariants 
which is particular efficient in treating two dimensional systems; in section 4, for sake of completeness, we 
recall systems admitting a second invariant which is a linear polynomial in the momenta; in section 5 we treat 
the case of the quadratic second invariant invariant; in section 6 we illustrate all known strongly-integrable 
solutions in the quadratic case and in section 7 we conclude. 



2 Hamiltonians with scalar and vector potentials 

We are interested in finding integrable examples of systems generated by a Hamiltonian function of the type 

U = -{pi + p 2 y ) + A x {x, y) Px + A 2 (x, y) p y + V(x, y), (1) 

where the function V is the ordinary "scalar" potential and A\ and A 2 are the component of a "vector" 
potential A in two dimension. Under the canonical transformation 

Px — ► Px + d x F, p y — > p y + d y F, (2) 

where F(x,y) is an arbitrary function, the Hamiltonian remains form invariant if A and V are changed 
according to 

A — > A + VF, (3) 
V — ► V + A- VF+±\VF\ 2 . (4) 
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However, the two quantities 

n{x,y) = k{d y A l -d x A 2 ) (5) 

and 

W(x,y)=V -i|A| 2 (6) 

are "gauge invariants" and can therefore be used to uniquely characterize the model system. The Hamiltonian 
to be worked on becomes then 

H = \iPx + A x {x, y)) 2 + l -{ Pv + A 2 (x, y)) 2 + W(x, y). (7) 

Q, the "curl" of the vector potential, has several physical interpretations: in astrophysical and celestial 
mechanical applications, it usually denotes an angular velocity field; it is a magnetic field in electrodynamics 
and plasma physics and so on. We remark that in general it is easier to attempt to solve directly for £1 and 
W. To recover the scalar V, one must have A x and A 2 and this can be another difficult problem [3]. 
We can first write the canonical equations provided by Q, 

x = p x + A u (8) 
y = P y + A 2 , (9) 
Px = -d x Ai p x - d x A 2 p y - d x V, (f 0) 

PlJ = -dyAl P X - 8yA 2 Py ~ 8yV, (ll) 

and then simplify them by exploiting the functions introduced above to get the equations of motion 

x-2fly = -d x W, (12) 
y + 2Qx = -d y W. (13) 

ft is readily verified that under the phase-space flow generated by (|12H13II , there exists a conserved function 
that is the first invariant of the system (Jacobi constant) 

J= ±(x 2 +y 2 ) + W. (14) 

In the investigation of the integrability properties of Hamiltonian , it turns out to be very helpful to 
work with complex variables. We perform then the canonical point transformation given by 

z = x + iy, p z =p= i(p x -ip y ), (15) 
z = x-iy, pg=p= h(j>x + iPy), (16) 

so that Hamiltonian Q turns out to be 

H = 2(p + $)(p+$) + W{z,z), (17) 

where the complex function 

§ = \{A x -iA 2 ) (18) 
has been introduced. In these variables, O is given by 

n = 23{9 2 $}, (19) 

where 3 denotes the imaginary part. Equations H15H16|I display a nice space-saving feature of using complex 
variables: even if an expression is not real, it is enough to write a single relation between complex functions 
(like, e.g. (|15fl ). The remaining information is provided by the corresponding complex conjugate expression, 
which we therefore do not write explicitly. For example, the canonical equations given by l|17|) are 

z = 2(p+*), (20) 
p = -2(p + ®)d z ®-2{$ + $)d x $-d z W, (21) 
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and the equations of motion corresponding to 1|12H13[1 are 



z + 2iilz = -2d g W. (22) 

The Jacobi constant now is 



J=±zz + W. (23) 



3 Polynomial invariants 



We are working with a Hamiltonian system with two degrees of freedom of which we already know an 
invariant, the Jacobi constant. In order to identify integrable systems in the usual Liouville- Arnold sense, 
we have to find a second independent phase-space function conserved along the flow. The standard direct 
method to solve the problem consists in making a suitable ansatz about this function and trying to solve 
the system of differential equations ensuing by the conservation condition. For several reasons, the ansatz of 
a polynomial in the generalized momenta is the most common 3 . It is well suited from the mathematical 
point of view, since it allows to get a system of PDEs in the coordinates only and is also well grounded on 
the basis of experience with already known integrable systems. 

Since we are looking for a real function, we make the following assumption 

M 

7 (M) =£(D fc (p + *)* + £ fc (p + *)*). (24) 

fe=0 

Although it is common, in the vector potential case, to see the invariant written in terms of the velocities, 
this is the correct interpretation of j(*0 as a phase-space function of the canonical variables. In order to 
satisfy the conservation of function (J22J, we impose that its Poisson bracket with the Hamiltonian vanishes, 

{/ (M) ,H} = (25) 

and try to solve for the complex functions Dk(z, z). In the presentation of the results, for easy comparisons 
with existing works, we will revert to the usual expressions in terms of the velocities, replacing the momenta 
in 12411 according to l|2U|l . In this case, conservation of the invariant can be checked by means of the condition 

^r = °' < 26 > 

along the solution of the equations of motion (1221 . 

Following the approach already used in the scalar case 01 0E], we consistently apply the trick of the 
energy constraint even in the present case. Here, with 'energy', we mean the Jacobi constant (|23|l . This 
method has the advantage of allowing the simultaneous treatment of "strong" invariants (the usual ones which 
are conserved for arbitrary values of the energy) and "weak" invariants (functions which are conserved only 
on some energy surfaces). In the papers cited above, we have moreover shown how the energy constraint 
simplify the structure of the system of PDEs that has to be solved. The essential remark is that, to identify 
the cases of strong integrability, it is sufficient that in the final results a subset can be isolated which is 
independent of the energy parameter, in the present situation the given value of the first invariant, let us say 
C. If we are interested in a "strongly-integrable system" , in the end we must get a solution independent of 
C. The procedure may appear more involute, but, at least in the scalar case, it reveals to be very effective. 

Formally, the idea of operating with the energy constraint is very simple. It consists first in introducing 
the "null" Hamiltonian 

n = n-c = o, (27) 

or, on the same footing, the "null" Jacobi invariant 

J = J -C =^zz-G = 0, (28) 
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where C as usual denotes the given fixed value of the Jacobi constant and 



G = C-W (29) 

is the so called "Jacobi potential" . Second, using in a consistent way the constraint 

2(p + $)(p+ $) = G, (30) 

wherever it appears in the computations. This essentially occurs when, implementing the conservation 
condition l|25|) , a polynomial relation in the generalised momenta appears and the constraint H3U|) is used to 
eliminate powers of the mixed variables (p + $)(p + $) in favour of powers of G. One minor shortcoming of 
this approach is that, in view of the explicit appearance of the function G, rather than simply the potential 
W, some of the coefficients of the invariant in general depend on the 'energy' parameter too: 

Dj = Dj{z,z;C), 0<j<M-2. (31) 

Therefore, to obtain the standard expression in terms of phase-space coordinates only, in the end we have 
to remember to perform the substitution 

C — W(p,p,M), (32) 

wherever the parameter C appears. In view of (|31|l . we see that this replacement does not affect the degree 
of the polynomial in the momenta. 

In practice, computing the Poisson bracket (|25[) . using the constraint, collecting the coefficients of the 
various powers of p + $ (they are accompanied by their complex conjugates) and imposing their vanishing, 
we get the system 

d- z D k - 1 +iknD k + ^d z (G k+1 D k+1 ) = 0, k = 0,l,...,M, (33) 

where it is implicitly assumed that Dj = for j < and for j > M. The set of equation l|33|l must be 
supplemented by the closure equations 

d- z D M = (34) 

and 

3t{d z (GD 1 )} = 0, (35) 

where -ft denotes the real part. For sake of space, we do not write the expressions of the standard direct 
approach in real coordinates and without the energy constraint. To compare with, we refer to Section 4 of 
Hietarinta |3| and recall the work of Hall where the study of weak invariants was first addressed and 
of Sarlet et al. |15| . where some wrong deductions contained in Hall's work were corrected. A systematic 
analysis of the cases with M = 1 (linear invariant) and M = 2 (quadratic invariant) was started in Ref.jZj 
and recently taken up again in Ref.|S]. A more general version of the problem concerned with quadratic 
invariants is mentioned in a different context in Ref.|ll)j. 

The first result one easily get with this approach is that equation (|34|l is readily solved as 

D M = D M (z), (36) 

that is, the leading order coefficient in the invariant is an arbitrary analytical function. This result, already 
known in the inertial case, still holds here. In Ref.jl] it was shown how, in the purely scalar case, the strong 
conservation condition restricts the form of this function. In what follows we will get analogous results when 
also a vector potential is present. 

One fundamental difference with what happens in the scalar case, is that equations for coefficients with 
even and odd indexes do not decouple. This fact is due to the Hamiltonian not be reversible in the present 
instance. System (|33f) is therefore very awkward to solve in the general case. In the next two sections we 
present the solutions in the linear and quadratic cases at fixed and arbitrary values of the Jacobi constant. 



5 



4 Linear invariants 



We start the investigation looking for systems admitting a second invariant which is a linear function in the 
momenta. The ansatz is 

= S(p + $) + S(p + $) + K, (37) 

where for the three coefficients we have used a notation which conforms with that in previous works. The 
system of equations ensuing from the conservation condition is the following: 

S- z = 0, (38) 
Kg + iflS = 0, (39) 
5R{(G5) 2 } = 0. (40) 

In order to compactify formulas, from hereinafter with the subscript we denote the partial derivative with 
respect to the corresponding variable. 

4.1 The general solution for linear invariants 

Equation (|38() agrees with l|34|l . confirming that S can be an arbitrary analytic function, 

S = S(z). (41) 

To complete the treatment, given an arbitrary S(z), we have to solve ()39|l and (|40|l . This task is more 
efficiently achieved by performing a coordinate transformation that trivialize the differential equations. Let 
us consider a conformal transformation z = F(w) to the new complex variable w = X + iY given by 

^ = F'(w) ee S(z(w)). (42) 
aw 



The explicit form of the transformation is then 



/dz 
W> ,43) 

and we have the relation between the differential operators 

A-=F>(w)±=S± (44) 
aw dz dz 

Multiplying (|40|l by the real factor SS, the content of the curly brackets can be modified in the following 
way: 

SS{SG) Z = S(SSG) Z = (SSG) W . (45) 
Introducing the "conformal" potential 

G=\F'\ 2 G=\S\ 2 G = SSG, (46) 

ca. (|40|l reduces to 

^{G w } = 0, (47) 

which is readily solved in 

G = g(Y), (48) 

where g is an arbitrary real function and, according to the definition of the coordinate transformation, Y is 
the imaginary part of w. 



G 



In the new variables w, eq.JSHJ), together with its complex conjugate, can be rewritten as 

K w = iti, (49) 

= -ih, (50) 

where the conformal field 

n = \s\ 2 n (51) 

has been introduced. The integrability condition for the real function K is 

St{K w } = 0, (52) 

with solution 

K = k(Y), (53) 
where k is another arbitrary real function. The conformal vector potential is then given by 

fi (54) 

An inversion of the coordinate transformation allows to express the solution in the original variables. The 
second invariant can be expressed as 

= di{S}x + 3{S*}y + K. (55) 
4.2 Linear invariants at arbitrary energy 

Equation (|40|) . in view of (|41|l and recalling the definition of G in Q29J1. can be rewritten as 

$t{S'(C-W)-SW z } = 0. (56) 

If we are interested in strong integrability, namely in an invariant which is conserved for arbitrary values of 
the Jacobi constant, equation (|56[1 must be independent of C. Therefore, it decouples in two independent 
equations: the first is 

5R{S"(z)} = 0. (57) 

The second turns out to be 

R{(SW) Z } = 0. (58) 

Ea. (|57Jl means that at arbitrary energy, we are no longer free in the choice of the coordinate transformation: 
we have to comply with this condition which actually imposes very strong limitations. It can be integrated 
to give 

S(z) = ikz + a, (59) 

where k is a real constant and a a complex constant. We can prove that it essentially allows only two kind 
of new coordinates: a) polar coordinates; b) rotated Cartesian coordinates. To show this, we first observe 
that we can exploit translations and scaling of the complex plane to further reduce the freedom contained 
in l|59l) . If k is not zero, a translation allows to put a = 0. A scaling allows then to put k = 1. We have 
then the two possibilities: 

a) S(z) = iz, F{w)=e iw , x = e~ Y cosX, y = e~ y sinX, (60) 

b) S(z) = a, F(w) =aw, x = aX - bY, y = bX + aY, (61) 

Case a) can be recognized as the transformation to polar coordinates. In fact, with the usual notation, they 
are defined as 

r = e _y , 8 = X. (62) 



7 



Solutions (|48H54|) . in view of (J2SJ, are 

W = W(r) = ^, K = k{r), ft = ^ (63) 

and the second invariant turns out to be 

1^ =i(zz-zz) + K = r 2 8 + k(r). (64) 

The problem is rotationally symmetric and the corresponding invariant is a generalisation of the angular 
momentum. 

Case b) can be recognized as the transformation to rotated Cartesian coordinates. Solutions H48I - I54JI ■ in 
view of then gives 

W = W(ay - bx), K = k{ay~bx) 1 r> = - fc ' (a ^~ bx ^ , (65) 
The second invariant now is 

= z + z + K = ax + by + k(ay ~ bx). (66) 

The problem is invariant under translation along the family of straight lines ax + by = const. These two 
solutions are already well known |7j, and the above procedure can be appreciated in its effectiveness. 



4.3 Examples of weakly integrable systems with linear invariants 

We may provide two interesting classes of weakly integrable systems admitting linear invariants. 

The first is obtained by the simple observation that, if we chose the level surface C = 0, it is no more 
necessary that condition Ij57(l be satisfied. Any analytic function S — S(z) provides a solution through the 
corresponding conformal transformation. If Y, as above, denotes the new coordinate 



S(z) 



then the solution is given by 



W= 9 -^», K = k(Y( X ,y)), = -^^l, (68) 

with g and fc arbitrary real functions. 

The second class of weakly integrable systems is obtained with the following trick. Let us consider the 
analytic function f(z) and consider then the conformal transformation l|42|l with S(z) given by 

' (69) 



c + /'(*)' 

with c constant. Recalling definition (|46|) . let us consider the "flat" conformal potential G = 1. In this case, 
relation l(4^|) . using ||BSJ|, gives 

G = C-W=j^ = c 2 + c(f + f') + \ff. (70) 

We can therefore interpret c 2 as the fixed value of the Jacobi constant 

c 2 = C (71) 
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and get as a consequence a C-dependent potential 

W(z,z;C)=-VC(f' + f')-\ff. (72) 

To complete the solution, we have to write explicitly the coordinate transformation generated by l|69[). that 
is 

■^r=cz + f(z). (73) 
Again, an arbitrary function k(Y), with 

Y = 3 {c* + /(*)}, (74) 
will do the work. For further details on these systems we refer to Ref.[TT]. 

5 Quadratic invariants 

We now look for systems admitting a second invariant which is a quadratic function in the momenta. The 
ansatz is 

I {2) = S( P + $) 2 + S(p + #) 2 + R(p + $) + R(p + $) + K, (75) 

where, besides S and the real function K, we now have to determine the complex function R. The system 
of equations ensuing from the conservation condition is the following: 

Sg = 0, (76) 

R s + 2iQS = 0, (77) 

K s + SG z + ±S'G + iMi = 0, (78) 

n{(RG) z } = 0. (79) 

where, in (|78J) we have already exploited l|76(l . that, as usual, embodies the fact that S is an arbitrary analytic 
function. 

5.1 Towards a general solution for quadratic invariants 

System l|76H79|l is much more difficult to solve than the previous linear case. Indeed we lack of a general 
solution. The main reason for this difficulty is that the coupling between G, R and f2 produces an integrability 
condition for K, through eq.J7S|) and its complex conjugate, that is a nonlinear PDE. However, we can 
implement the strategy to arrive as close as possible to a general solution and, what is of great importance, 
we can solve the problem in the strongly intcgrablc case, developing an effective way to construct solutions. 
This time we use a conformal transformation to the complex variable w = X + iY given by 

— = F'(w) ee V / 5(^H) (80) 



so that the explicit form of the transformation is 



dz 



(81) 



Introducing the conformal potential 

G =\F'\ 2 G = \S\G = yfs^G, (82) 



and the new complex function 



~ R = Ts- (83) 
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ea. (|79|l keeps its form in the transformed coordinates 

$i{(RG) w } = 0. (84) 

The solution of this equation is 

RG = HC W , (85) 

where K is an arbitrary real function. 

Let us now multiply both sides of l|77|) by 

\/§. TO 

Using l(83j) and introducing the conformal field 

n=|£|n, (87) 

cq. (17711 transforms into 

R a + 2tQ = 0. (88) 
Since the conformal field is real, the solution of l|88|) is 

R = -U£ w , (89) 

where £ is another arbitrary real function. The factor —4 appears for later convenience. In this way, the 
conformal field is given by 

n = 2£ wiB . (90) 

At this point, to get the general solution of the unified treatment of weak and strong integrability, one 
should try to solve the integrability condition for ea. l78|) in the light of the results l|85H89[) . Let us write the 
integrability condition for K which, computing K zz = K zz from J7SJ, is 

3{5"G + 3S'G Z + 2SG ZZ + 2i(QR) z } = 0. (91) 

In the transformed coordinates, this becomes 

${G ww }+$l{{Rh) w } = 0. (92) 

^From (|S5|> and (|55|l we have 

i?0 = ^RRu, = j{R 2 )w — — KG" 2 ^^ 2 )^ (93) 

so that (|§2"|l can be rewritten as 

%{G WW } = \^{i{G- 2 1C^ w } (94) 

We can try to solve this equation after specifying JC. Unfortunately, the equation is highly nonlinear except 
in the rather trivial case when K, is constant, so it appears to be very difficult to solve it in the general case. 
We direct our attention attempting to solve the more limited but fundamental case of strong integrability. 

5.2 Quadratic invariants at arbitrary energy 

Starting again with l|79|l . we see that, in view of definition 129(1 . it can be rewritten as 

$t{R z (C-W) -RW Z } = 0. (95) 

If R is independent of C, as is the case under study, in order to satisfy this equation at arbitrary values of 
the energy, it decouples in two independent equations: the first is 

n{R z } = 0; (96) 
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the second 

$t{RW z } = 0. (97) 
Eq. (|96[) is solved by introducing the arbitrary real function 77 so that 

R = -4irjz. (98) 

This allows to solve equation l|97|) for W in the form 

W = W{rj), (99) 

that is W is, at this stage, an arbitrary function of the argument. 

Comparing (|9"5|l with (|S9")l and taking into account definition {H3J|, we get the relation 



SS£ W = \S\£ W = 7] W . (100) 

Solving the integrability condition for 77 with a given form of 5* allows to find £ (and viceversa): equation 
(|100fl is deceivingly simple; we will see later that it is of some concern. Here we remark that the route 
followed to treat ea. H79l) reverted to the original physical coordinates z,z, since they lead to the simple 
result in l|98|) . However, ea. (|90|l above is valid in general: it is the simplest way in which we can solve for 
the vector potential and is expressed in the new coordinates w,w. This result and the development below 
show how working in the new coordinates is advantageous in this context too. 

In the meantime we have to determine the forms of S imposed by strong integrability. Examining the 
integrability condition (|91|l . the usual condition of independence of the results from C imposes the following 
constraints on the form of S 

%{S"(z)} = 0. (101) 

This result, which is valid in the scalar case also, is a natural extension of what found in the linear case. 
Coming back to ea. (jlOOfl . we note that, using real coordinates, we have 

\S\ix = Vx, (102) 
= ~VY- (103) 

It can be proven ^1] that condition (|101Jl is equivalent to 

\S\xy = 0, (104) 

so that we can also write 

\S\=A(X)+B(Y), (105) 

with A and B determined by the specific form of S. Therefore, eas. l|lU2(llU3fl generates the following 
differential equations for the functions rj and £ 

A'r) Y + B'r)x-2(A + B)r)xY = 0, (106) 
A'Zy+B'Sx+2(A + B)Zxy = 0. (107) 

In view of (|98|l . that can be rewritten as 

|5|E = -4^ (108) 
and of (|87|l . integrability condition l|92() can be written in the form 

Z{G ww } + \S\n{m w }=0, (109) 

or, in real coordinates, 

GxY + 2( Vx Q Y -r lY n x ) = 0. (110) 
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If we want to exploit the result established by l|99|) . we have to resume the physical potential W through 

G=\S\(C-W). (Ill) 

Using l|105|) we have 

Gxy = -A'Wy - B'W X ~(A + B)W XY (112) 

and so [jlltty is 

(A + B)Wxy + A'Wy + B'W X = 2(r?x0y - m^x)- (113) 

But lO tells us that 

W x = nxW'{r,), W Y = r, Y W'{r,), W XY = VxVyW" + vxyW, (114) 
so that (|113fl becomes 

W{n) + 3^-W{ v ) = f- - — ) ■ (us) 

rj X r/Y A + B \rj Y r)x ) 

This is the best form we attain to express the integrability condition for the potential: we see that it implies 
that 

J^=$(77), (116) 
VxVy 

where $ is arbitrary. 

Therefore, the strategy to find strongly integrable systems with scalar and vector potentials supporting 
quadratic invariants is: to chose a suitable S in the class determined by condition H1U1|) (and therefore A 
and B) and solve (|i06lli07|> to find Tf(X,Y) and £(X,Y); to solve © to find 

^ = \(^xx+^yy) (117) 



and use 1)87(1 to find 

and, finally, try to solve l|115|) for 
taking into account nitty . 



q{x ' y) = aTb (118) 



W = W (r)(X,Y)) (119) 



5.3 General form of the quadratic invariant 

In the next section we apply the strategy delineated above illustrating how already known and new integrable 
systems are determined. We end this section with a closer look at the structure of the quadratic invariant. 

In the new variables, w,w, together with the conformal transformation, it is natural to introduce the 
new time variable r such that 

dT = Wf (120) 

We can use the apex to also denote the derivative with respect to r without risk of confusion with other 
derivatives with respect to coordinates. Equations of motion ll'l'li then assume the form: 

w" + 2m W ' = 2G W . (121) 

In view of above positions, the quadratic invariant is most simply expressed as 

I (2) =ft^(w') 2 + Rw'^+K, (122) 
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or, in terms of real variables, 



/(2) = I 
2 



{X'f~(Y'f - 2{UY' + £ Y X') + K. 



(123) 



We recall that the scalar K, to be found by integrating the transformed version of J7SJ), namely 

Ku, + G w + iQR = 0, (124) 

is a function of the form 

K = K(w, w; C). (125) 

In the present case of strong integrability, in view of (|1U5[) . it turns out that this function can be expressed 
as 

K = k + C(B — A), 
and k can be found by integrating the system 

k x -Wx+m x = 0, 
k Y + W Y - 4f2£y = 0, 

where W = {A + B)W . From using l(T2^|) and 



l (X') 2 + (y') 2 
2 A + B 



W, 



the general form of the invariant is then 

1 



j(2) 



B {X' f - A (Y'f - 2(CxY' + &X 1 ) + (B — A)W + k. 



A + B 

Observing that the relations between velocities in the two gauges are given by 

X' = ft{F'}x + %{F'}y, 
Y' = -%{F'}x + n{F'}y, 



(126) 

(127) 
(128) 

(129) 
(130) 



(131) 
(132) 



we can eventually transform the invariant in Cartesian coordinates. Transformation rules automatically 
account also for the change of the time variable according to 112U|) . 



6 Solutions with quadratic invariants 

Let us recall the coordinate systems given by the condition for the existence of strong quadratic invariants. 
We observe that the function S we have to use is obtained by integrating I|1U1|) so that 



S(z) = cz 2 + f3z + a, 



(133) 



where c is a real constant and /3, a complex constants. Exploiting the freedom of making translations, 
rotations and scaling in the complex plane, we have the following four inequivalent cases (for further details 
we refer to Refs.gl EI]): 



a) S(z) 

b) S(z) 

c) S(z) 

d) S(z) 



x = aX -bY, y = bX + aY, 



a, F(w) = aw, 
z 2 , F(w) = e w , 
4z, F(w)=w 2 , 

z 2 + A 2 , F(w) — Asinhw, x = A sinh X cos Y, y = A cosh X sin Y. 



x = e x cos Y, y = e x sin Y, 
x = X 2 - Y 2 , y = 2XY, 



(134) 
(135) 
(136) 
(137) 
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Case a) gives the rotated Cartesian coordinates (standard Cartesian coordinates if a = 1). Case b) gives 
again the polar coordinates. Case c) gives the parabolic coordinates which are also referred to as Levi-Civita 
coordinates (the factor 4 in the definition of S appears just for this reason). Finally, case d) produces the 
elliptical coordinates. For sake of completeness, we list, for each of the four cases, the conformal factor and 
functions A and B introduced in (|105|l : 

a) 15*1 = 1, A = S=i, (138) 

b) \S\ = e 2X , A = e 2X , B = 0, (139) 

c) \S\ = 4(X 2 + Y 2 ), A = 4X 2 , B = 4Y 2 , (140) 

d) \S\ = A 2 (sinh 2 X + cos 2 r), A = A 2 sinh 2 X, B = A 2 cos 2 Y. (141) 

6.1 Cartesian coordinates 

Let us start with the simplest case a), the Cartesian one. As the structure of the solution will show in the 
end, there is no need to work with rotated coordinates, therefore we can put a = 1 and make the trivial 
identification X — x and Y = y. Eas. H102H103|l are 

Vx = (142) 
% = (143) 



The solution is 



C = f(x)+g(y), (144) 
V = f(x)-g(y), (145) 



with / and g arbitrary. From l|117|) and l|118|l . the vector field is 

n = \(f"+9") (146) 

and the R function 

R=-2(g' + if). (147) 

Eq.dl is 

W xy - (g'f" + f'g'") = 0. (148) 
^From (inn we have that W = W(f(x) - g(y)), so that 

W xy = -f'g'W" (149) 

and i|148|) becomes 

fin in 

W" + —— + —— = 0. (150) 
/' 3 

A further differentiation by x and y gives 

1 d f" 1 d a'" 

^r-?lr fc - r (151 » 

with a real constant. We then get two equations for / and g 

f" = af 2 + bf + c, (152) 
g" = -ag 2 + dg + e, (153) 
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in agreement with the work of Dorizzi et al. [7], that can be integrated for various choices of the constants 
a, 6, c, d, e. The potentials in terms of / and g are 

W = ^g-ff- b -±l( g -ff + ic -e + m)(g-f), (154) 
2Q = a(f - g 2 ) +bf + dg + c + e, (155) 
where m is another constant. Integrating eas. (|127H128|l gives then 

K = afg(f + g) + (b - d)fg - a(f + 9 3 )-\ (36 + d)f + \ (3d + b)g 2 - (3c + m + e)f + {c-m + 3e)g. (156) 
The fundamental example of the harmonic oscillator is obtained with the choice 

a = b = d = 0, (157) 



so that 



/ = \cx 2 7 (158) 
9 = \ey 2 (159) 



and the potentials are 



W = \{c- e + m){ey 2 - cx 2 ), (160) 
20 = c + e. (161) 

We remark that, in this case, the vector field is a constant angular velocity and that the isotropic oscillator 
only exists when the angular velocity vanishes. The second invariant in this case is 

/< 2 ) = l(i 2 - y 2 ) - 2(cxy + eyx) + ^[(c-m + 3e)ey 2 - (3c + m + e)cx 2 ] . (162) 
6.2 Polar coordinates 

In case b), the polar coordinates, we recall that now 

r = c x , 9 = Y (163) 
and remark the difference of this relation with respect to (|62|l obtained in the linear case. We have 

\S\= A = c 2X =r 2 , B = 0. (164) 

Ea. p()6(l becomes 

2e 2X {n XY ~ VY ) =0, (165) 

which is solved by 

7 ? = c x g(Y) + f(X), (166) 
with / and g arbitrary. From i|l(JU[l we then have 

£ = -e x g+ fe- 2X f'(X)dX. (167) 

/From (|117[) and !|118|) . the velocity field is 

2n = e~ 2X A£ = e-«(/" - 2/') - <r* x {g" + g). (168) 
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Condition 1)11611 says that 

Vx = c x g(Y) + f\X) = $(e x g(Y) + f(X)). (169) 

There are two ways in which l|169f) can be accomplished: the first is g = const and it is easy to see that in 
this way we are actually taken back to the spherically symmetric case to which actually a linear invariant is 
associated. We can take r/(X) and W(X) arbitrarily whereas, to find the vector field, we observe that 

£'P0 = e~ 2 VPO (170) 

so that 

2fl = c- 2X e = G ~ AX {ri" - 2r/) = e~ 2X (e~ 2X j)' . (171) 

To show that this solution is equivalent to that obtained in the linear case, we compute the scalar k by 
integrating (|127fl to get 

k^(X) = c 2X W-(e- 2X V ') 2 . (172) 
Observing that, in polar coordinates, relation (|120|l between old and new time is 

dr = % (173) 

and using 

Y' = r 2 (), (174) 

expression ()130JI becomes 

j(2) = _ r 4^2 _ 2 jQ _ (1?5) 

In the linear case, recalling solution l|64(l . we got: 

= r 2 9 + k [1 \r), (176) 

where, from 163fl . 

k^\r) = 2 J rfldr = 2 J e 2X ndX. (177) 

Using 1)171(1 this becomes 

fc(i)(r) = e - 2X r/ = r-V, (178) 
so that, comparing l|175|) with l|176|) . we see that 

=-(/W) 2 . (179) 

The second possibility of satisfying condition (|169|l is to have / = const and the constant, without any 
loss of generality, can be put equal to zero. Therefore, we have 

r, = c x g(Y) (180) 

and 

£ = -e~ X 9- (181) 
To find the potential, we have to solve l|115|l which, using (|180H and (|168fl . is 

W"(v) + -W'( V ) = -- J-^ (g"'g + 3g"g' + 4g'g) . (182) 

The simplest way in which the right hand side too is a function of ry is given by the condition 

g"'g + 3g"g' + Ag'g = a ^, (183) 
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with a constant, so that (|182[) becomes 



W"(ji) + -W'(ri) + 4 = (184) 

with solution 

"-?-■§? < 185 > 

To find explicit solution we have to determine the possible forms of g which can be got by integrating (|183[1 
twice 

(5') 2 +ff 2 + A + ^=^ ( 186 ) 
4 5 4 g A 

with 7 and 5 constants. All this is in agreement with the results in McSween and Winternitz 8 to which 
we refer. 

Since in the approach with the energy constraint, the coefficients appearing in the invariant must be 
supplemented by a further substitution in order to get the final physical expression, for sake of completeness, 
we work out a specific example. To ease the comparison with the results in Ref.[5], we use explicitly polar 
coordinates as in (|163fl . The simplest case in the solution above is given by the choice a = 0. Ea. (|186l) has 
solution 

g{9) = Va + bcos26 (187) 

with the constants a and b given by 

7 = a 2 -6 2 , 5 = 2a. (188) 
^From Ijl85|l and Hlfi8|l . the potential and velocity field are given respectively by 

P 

W= r 2 (a + bcos29) (189) 

and 

Q = 1,2 - a " (190) 

2r 3 (a + 6cos2(9) 3 / 2 ' V ' 

To find the invariant, we can use l|13Ufl . Eas. H27H128fl can easily be integrated to get 

k = -r 2 W — r~ 2 g(g + g"). (191) 



so that expression (|130fl becomes 



1^ = + 2(rg8 - rg 1 ) + % + g(g+ / . (192) 

In the specific example of l|187(l we get 

/(2 ) = r 4^2 2bs[n9 r + 2Va + bcos26r6 + ± ( + 2(3] . (193) 

Va + 6cos26» a + &cos26»V r 2 ' J y ' 

We remark that the small discrepancy in (|193J) with respect to the analogous expression reported in Ref. [B] 
is due to a difference of a factor 2 in the definition of O as can be seen from ea. ()19()(l . 

6.3 Parabolic coordinates 

We use the Levi-Civita representation of the parabolic coordinates introduced in case c) above and mostly 
used in celestial mechanical applications. As a first step, the general strategy we depicted at the end of 
subsection 5.2 prescribes to determine the functions ij and £ satisfying equations HIUGHIOT]) . According to 
(|14(jp . we have 

A = AX 2 , B = AY 2 . (194) 
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so that ea. (|106|) becomes 

X m + Yi lx -{X 2 + Y 2 )r, XY = Q. (195) 
A fairly general solution of this PDE can be represented in the following form 

V(X,Y) = [ F(a)y/{X 2 + a)(Y 2 ~a)da + [ G(a) i^l±I^T 



—da, (196) 

J J {(X 2 + a){Y 2 - a)f /2 

where a is an arbitrary real parameter and F and G two arbitrary smooth functions. 

The subsequent steps should consist in determining the functions £ and fl and finally to solve for the 
potential W(rf). However, we observe that it is not easy to satisfy the constraint (|116|) with a too general 
expression for rj: therefore, we first find a suitable form of this function and then proceed as above. A simple 
but non trivial possibility is that given by the position 

F = cS(a-b), G = (197) 

in (|196fl . where 5(a) is the Dirac function. This choice gives the simplest separable solution of the differential 
equation l|195|) . Therefore, we have 



V(X, Y) = c^(X 2 + b)(Y 2 -b), (198) 
which is easily shown to satisfy condition i|116|l in the form 

Vxy 1 



rixm V 

With solution (|198f) . equations (|10(jHl(J , /fl are readily solved for £ 



(199) 



^(X,y) = 7 arctan A /^— (200) 



I X 2 + b 

4 V V 2 



so that, from l|117|) . the conformal vector field is 

((X 2 + b)(Y 2 - b)) 3/2 



~ be X 2 + Y 2 

( -> - ... . . (201) 



^From this, ea. (|118fl gives 



This result suggests to put 



n = ? = ^i. (202) 

4(A 2 + Y 2 ) 32 T? 3 { ' 



c = 2 (203) 
in order to simplify formulas and we have the physical vector field from 

- A V (X, Y) = 2^/{X 2 + b)(Y 2 -b). (204) 

An important consequence of this result is that, together with W, £1 too depends on the coordinates only 

through rj. This implies that the right hand of ea. (|115|l for W vanishes, so that, taking into account 1199|l . 
we get simply 

W"(r)) + -W'(r))=0. (205) 
V 

The solution for the scalar potential is then 

W = 4 (206) 
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and we are actually led to a situation analogous to that encountered above in the example examined in polar 
coordinates. It is interesting to remark that potential l|206|l is separable if considered in the purely scalar 
situation, since, using l|198|) . we get 



W 







1 



4{X 2 + Y 2 ) \X 2 + b 
Reverting to Cartesian coordinates through 



Y 2 -b 



X = 



r + x 



Y = 



the scalar and vector potentials are respectively given by 

/3 



W 



and 



CI = 



y 2 - 4b(x + b) 
b 



2(y 2 -4b(x + b)) 3/2 ' 
The second invariant using parabolic coordinates is 

j(2) 1 



, Y 2 (X'Y - X 2 (Y'Y 
X 2 + Y 2 L y ' y ' 

Y(X 2 + b)X' - X(Y 2 - b)Y' 



(X 2 + Y 2 )^(X 2 + b){Y 2 -b) 
b + 8[3(Y 2 -X 2 -b) 
4(X 2 + b)(Y 2 ~ b) ' 



(207) 

(208) 

(209) 
(210) 



(211) 



where the r time variable is used and, therefore, it is conserved along the solution of the transformed 
equations of motion of the form l|121|) . Using relations I|131II132'|) between velocities in the two gauges, 



X' 

V 



y/2(r + x)x + y/2(r-x)y, 
y/2(r + x)y - y / 2(r-x)x, 



we can transform the invariant into Cartesian coordinates 

yx + 2by 



j( 2 ) = ( yx - xy)y + 



(x + b) 



2^/y 2 - 4b(x + b) 4(y 2 -Ab{x + b))' 



(212) 
(213) 



(214) 



6.4 Elliptical coordinates 

It turns out that even in elliptical coordinates it is possible to find non trivial solutions with a structure 
closely related to that seen in the examples detailed above in polar and parabolic coordinates. In addition, 
in elliptical coordinates it also exists a solution with a constant vector field (constant 'angular velocity'), 
whereas this possibility appears to be absent in parabolic coordinates. 

A common feature of all these cases is that both the scalar and vector potentials, W and fl, can be 
expressed only in terms of rj. This implies that the right hand of cci. (|115fl for W vanishes, so that, taking 
into account l|ll(j|) . this common feature is embodied in the equation 

W"{rf) + §[ri)W'{ri) = 0, (215) 

with suitable $. Analogously to what seen in the polar cases, other solutions may very well exist but they 
are not so easy to find. 
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According to (|141f> . we now have 

A = A 2 sinh 2 A, B = A 2 cos 2 Y. (216) 

In this case ea. H106|) becomes 

sinh X cosh Xr\y — sin Y cos Yj] X - (sinh 2 X + cos 2 Y)t] X y = 0. (217) 
A simple solution of this PDE is the following 

rj(X, Y) = aA 2 (sinh 2 X - cos 2 Y) + feA 4 (sinh 2 X + cos 2 Y) 2 , (218) 
with a and b constants. Using this solution, equations H106l - I107|) give for £ 

£{X,Y) = a In [A 2 (sinh 2 A + cos 2 Y~)] + i&A 2 (cosh 2X - cos 2Y). (219) 

Computing the conformal vector potential from (| 1 1 T|> - we see that the first solution is trivial since its 
Laplacian vanishes giving a null vector field. We remark that the analogous phenomenon also occurs in 
the parabolic coordinates case, where it is possible to find several additional solutions of equation (|195fl 
generating a vanishing vector field. 

The second solution appearing in l|219l) instead gives 

Q = bA 2 (cosh 2 A + cos 2Y) = 2b A 2 (sinh 2 A + cos 2 Y). (220) 

We get that the conformal field is proportional to the conformal factor 15*1 so that, from l|118|) . we have that 
the vector field is constant 

= 2b=VL . (221) 
Condition (|llfi|> for the non trivial solution Ij218|) with a = is 



Vxy 1 



TjxrjY 2?7 

so that ea. (|115|) for W is 



(222) 



2 V 

whose solution is 



W"(v) + irW'to) = 0, (223) 



W = ——, (224) 

or, using again 1218|l . 

W(X, Y) = ^= ^ — , (225) 

1 6 1 A 2 (sinh A + cos 2 F) 

where a and f3 are arbitrary constants. This solution is the simplest separable scalar potential in elliptical 
coordinates, which is therefore shown to be integrable also in a uniformly rotating system: this result is 
already known for its application in the modelling of rotating galaxies |16| . The second invariant 113UI) using 
elliptical coordinates is 

r (2) - [cos 2 F(A') 2 - sinh 2 (r') 2 ] + ^(sin 2FA' - sinh2Xr') 



P Z , - ^o(sinh 2 A + cos 2 Yf 



sinh 2 A + cos 2 Y 1 v " ' v " ' J ' 2" 
sinh 2 A — cos 2 Y 
sinh 2 A + cos 2 Y 
Reverting to Cartesian coordinates through 

2A 2 sinh 2 A = r 2 - A 2 + ^(r 2 + A 2 ) 2 - 4A 2 ?/ 2 , (226) 

2A 2 sin 2 Y = r 2 + A 2 - ^(r 2 + A 2 ) 2 - 4A 2 y 2 , (227) 
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the scalar potential becomes 

W(x, y) = 13 (228) 

^r 4 + 2A 2 (x 2 -y 2 ) + A 4 

and, using relations (|131H132*)) between velocities in the two gauges, we can transform the invariant into 
Cartesian coordinates 

1^ ={yx - xyf - A 2 x 2 + (A 2 - r 2 ) (2fl yx - W(x, y)) + 2fl (A 2 + r 2 )xy+ 

tt 2 (r 4 + 2A V - y 2 ) + A 4 ). ( } 

We pass now to investigate a more complex class of systems. A general solution of equation l|217(l is 

r](X,Y)= J F(a) sj (A 2 sinh 2 X + a) (a - A 2 cos 2 Y)da+ 

(A 2 sinh 2 X + A 2 cos 2 Y) 2 , ( 230 ) 

G(a) ™ da, 

((A 2 sinh 2 X + a)(a - A 2 cos 2 Y)) 

where a is an arbitrary real parameter and F and G two arbitrary smooth functions. Comparing this solution 
with that of l|196(l we can guess analogous developements. Therefore we try with the simple solution 

F = 8(a-b), G = 0, (231) 

which corresponds to the simplest separable solution of the differential equation (|217|l . Therefore, we have 



r]{X,Y) = y(A 2 sinh 2 X + 6)(6- A 2 cos 2 y), (232) 

which satisfies condition (|116H in the same form as in (|199|) . With solution l|232() . equations H106H107|) give 
the following expression for £ 



so that, from l|117|) . the conformal vector field is 

bA 2 {b- 

((A 2 sinh 2 X + b)(b - A 2 cos 2 Y)) 



~ = &A 2 (6-A 2 )(sinh 2 X + cos 2 y) m) 



^From this, ea. (|118|l gives 

Q _b(b-A 2 ) 

A 2 (sinh 2 X + cos 2 r) ~" 2rf 

As remarked above, from this results it follows that the right hand of ea. (|115fl for W vanishes, so that, taking 
into account l|199(l . we get again l|205|) so that the solution for the scalar potential is 

0_ = P_ 

V 2 (A 2 sinh 2 X + b)(b — A 2 cos 2 Y) ' 



0= 5 = - 5-^-. (235) 



W = ^= ■ 7TT, , ,w, -o (236) 



Even potential l|23()l) is separable if considered in the purely scalar situation, since it can be written in the 
form 

W = t ( I I ^ (237) 

A 2 (sinh 2 X + cos 2 Y) \b - A 2 cos 2 Y A 2 sinh 2 X + bj 

Using the explicit coordinate transformation H226H227|l . the scalar and vector potentials are respectively 
given by 

W = by 2 + {b-A 2 ){x 2 + b) (238) 
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and 



a = Kb ~ A2) (239) 

2(by 2 + (b- A 2 )(a; 2 +6)) 3 / 2 ' V ' 



Finally, the second invariant using elliptical coordinates is 

cos 2 ypr) 2 -sinh 2 (y') 2 



j(2) 



sinh 2 X + cos 2 Y 



, „ Tr b + A 2 sinh X . . , „ v / b - A 2 cos 2 Y . 

+ 5 sin2rW- — — X'-sinh2XW 5 — Y 

sinh 2 X + cos 2 Y \ \ b-A 2 cos 2 F y b + A 2 sinh 2 X 

b(A 2 -b) + 2/3(A 2 (cos 2 F - sinh 2 X) - b) 
+ (A 2 sinh 2 X + b) (b - A 2 cos 2 Y) ' 

whereas using Cartesian coordinates is given by 

TO - (,* - *„? - AV + 2 W ^' - *f + K» - ^Wy^±a . (240) 

^foy 2 + (6 - A 2 )(x 2 + 6) ojr + (0 - A^)(a;^ + 6) 

7 Concluding remarks 

We have investigated Hamiltonian systems with vector potentials admitting a second invariant which is linear 
or quadratic in the momenta. In our approach, weak and strong invariants are treated in a unified setting 
where the strong invariants emerge as special cases. As for scalar potentials, the integrable systems can be 
greatly simplified by introducing certain standardized coordinates, as given in (|134I137|I . It is a striking 
result that these standardized coordinate systems for systems with strong invariants exactly coincide with 
the classical separable coordinates for scalar potentials. 

This work is an extension and improvement of the approach to integrable vector potential Hamiltonians 
which was proposed in jllj . However, there still remain issues which need clarification. In particular, it 
should be possible to obtain a better understanding of the integrability conditions, especially the role of the 
condition (|116fl for the structure of the strongly invariant case. 
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